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ABSTRACT 


Yen, Wen-Yo. Ph.D., Purdue University, December 1977. 
Effects of Dynamic Aeroelasticity on Handling Qualities 
and Pilot Rating. Major Professor: Robert L. Swaim. 

Pilot performance parameters, such as pilot ratings, 
tracking errors, and pilot comments were determined for a 
longitudinal pitch tracking task using a large, flexible 
bomber with parametric variations in the undamped natural 
frequencies of the two lowest frequency symmetric elastic 
modes. This pitch tracking task was programmed on a fixed 
base simulator with an electronic attitude-director display 
pitch command, pitch angle, and pitch error. The results 
of this study indicate that low-frequency structural 
flexibility can significantly affect the handling qualities 
and pilot ratings in the task evaluated. 
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CHAPTER 1 
INTRODUCTION 


1 . 1 Handling Qualities and Pilo t Rating 

The way an airplane in flight responds to control 
inputs by the pilot is referred to as the handling qualities 
of the vehicle. If the handling qualities ere not good, 
the pilot must devote more of his attention to flying the 
airplane and less to other mission activities, such as 
weapon delivery and air-to-air combat. In other words, 
the worse the airplane handles, the greater the pilot 
workload. In simple terms, handling qualities relate to 
the ease or difficulty in attempting to maneuver or control 
an aircraft. The study of handling qualities is 100% 
pilot oriented. 

Although all the individual factors contributing to 


handling qualities may fall in the satisfactory or accept- 
able ranges of their respective boundaries, their composite 
effect may still not produce a completely satisfactory 
airplane in some cases. Proof of the total system handling 


qualities is usually required by ground-based simulation 
and/or flight test of the actual airplane. A tool that 
is widely used in such evaluations is the Cooper-Harper 
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Pilot Opinion Rating Scale (reference 8). It is a ten- 
point numerical scale whereby the pilot can express his 
evaluation or opinion of how well the airplane flys or 
handles. 

1 . 2 Handling Qualities History 

There were only a few attempts to design an aircraft 
for good handling qualities before the second world war 
era. Before then the design of aircraft concentrated pri- 
marily on aircraft performance goals. Stability and con- 
trol characteristics were not well understood. It was 
known* however, that static and dynamic stability character 
istics and control system design determined the handling 
qualities - good or bad - of an aircraft. However, since 
that era much has been done to analyze the handling quali- 
ties of existing aircraft and to develop theories that 
allow the stability and control engineers to design an 
aircraft that handles well. 

Presently, designers mainly consider the rigid air- 
craft to estimate the handling qualities and pilot rating 
of the airplane. But large airplanes such as transport 
category airplanes and bombers need to be as light as 
possible due to many reasons which include fuel conser- 
vation. This means a sacrifice in structural rigidity. 
Airplane flexure causes additional aerodynamic loads which 

in turn cause additional flexure, etc. Also coupling 
occurs between the elastic modes and the rigid body motion 
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as the gyros sense the flexure motion and the rigid body 
motion. 

To point out the near complete lack of knowledge as 
to how dynamic aeroelastic modes affect airplane handling 
qualities, we quote the only reference to their effects 
contained in MIL-F-8785B, Military Specification - Flying 
Qualities of Piloted Airplanes (reference 2) . The four- 
line statement merely says: "Since aeroelasticity , control 
equipment, and structural dynamics may exert an important 
influence on the airplane flying qualities, such effects 
should not be overlooked in calculations or analyses dir- 
ected toward investigation of compliance with the require 
ments of this specification." The specification is con- 
cerned only with desirable ranges of values on rigid-bcdy 
dynamic response parameters. It seems quite possible that 
the handling qualities and pilot rating could be signifi- 
cantly affected by elastic modes, particularly in case 
some modes have low natural frequencies. It is not at all 
clear that the handling qualities should be specified by 
rigid-body dynamic parameters when such rigid and elastic 
mode interaction is present. The pilot will likely not be 
able to discern, for example, how much of a given pitch 

angle response to command input is due to rigid-body and 
how much is due to low frequency elastic modes. These 

interactions certainly affected the pilot's assessment of 
the handling qualities in the work reported herein. 
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Therefore, some study of these elastic effects and revised 
design standards must be inaugurated for large flexible 
vehicles. To date no design criteria exist for handling 
qualities in terms of control system specifications. It 
is hoped that this thesis will be a basic step in that 

direction. 

1.3 Objectives 

The longitudinal dynamics of the Rockwell B-l bomber 
aircraft was simulated on an analog computer. The two 
lowest frequency symmetric elastic modes were included in 
the math model. This was sufficient to get an apprecia- 
tion for the effects on pilot rating, where the elastic 
mode characteristics were varied between runs. 

A fixed base, pilot in the loop simulation was 
performed to examine the handling qualities and pilot 
ratings with variations in each of the undamped natural 
frequencies of the two elastic modes. Specifically, the 

following questions were addressed: 

1. What are the relative pilot ratings in a pitch 
tracking task as a function of lowering elastic 

mode frequencies? 

2. Which frequency combinations produce the poorest 
pilot performance? 

3. Which of these two elastic modes has the most 
significant effects on the handling qualities 
and pilot rating? 
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Pilot pitch tracking performance, Cooper-Harper 

rating, and pilot comments were used to answer these 
questions . 

4 Thesis Organization 

Chapter 2 mathematically establishes the task. The 
aircraft dynamic equations are derived and the simplifi- 
cations implicit in the perturbation technique and 
aeroelastic effects are discussed. In addition, a flight 
director equation and a pitch angle error equation, 
including their functions in this study, are presented. 

Chapter 3 discusses how and why the eight cases are 
chosen for this study. The dynamic characteristics of 
eight cases are presented. 

Chapter 4 discusses the experimental objectives, 
apparatus, and methodology for the fixed base, pilot in 
the loop simulation. 

Chapter 5 gives the results of this simulation study. 

The eight cases are compared and the difficulty of some 
cases is discussed. 

Chapter 6 concludes the thesis by reviewing the 
results, outlining areas that require additional study, 
and recommending further activities. 
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CHAPTER 2 

EQUATIONS OF MOTION 


2.1 funeral Des c ription of the Aircra ft 

The aircraft simulated in this study is the B-l 
bomber, an American built, large flexible aircraft 
(Figure 2.1). Rockwell international examined this 
aircraft which they designed (reference 3) and found 
they were able to fly a sea level zero degree trim pitch 
angle at 562.2 knots ,943 feet per second airspeed) with 
the elevators up 6 degrees. It is from these studies 
that much of the aircraft data (mass, stability denva- 

tives, etc.) was obtained. 

There is a reason the B-l bomber was used in our 

studies. This airplane is large and flexible enough to 

exhibit the desired effects between fuselage elastic 

motion and rigid-body motion. 

2 2 Airframe Equations 

In order to write equations of motion one must first 
establish a coordinate system. It is common practice to 
attach a body-fixed coordinate system to the aircraft 
and write equations with respect to it. By doing so a 
moments of inertia will remain constant. The stabi ity 


I 




axes convention is to let the X axis point forward, the 


Y axis point out the right wing, and the Z axis point 
downward. 


I. Rigid Body Equations 

Assuming constant mass: 



m ( 


dv 

3t 


+ w x V) 


B 


M 



= ^ + a) x H 

dt B 


( 2 .] 


( 2 .; 


where : 


dv 

dt 


I 


is the rate of change of the aircraft 


velocity vector as seen from inertia coordinates. 


I is the rate of change as seen from aircraft 

B 

body axes. 

ui is the angular velocity of the aircraft. 

H is the aircraft angular momentum vector. 

F is the vector sum of all external forces. 

M is the vector sum of all external moments, 

m is the aircraft mass. 
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These two vector equations can be written in terms 
of their components to produce the following six 


equations ; 

G x + F x » m(U + QW - RV) (2.3) 

G y + F y « m(V ♦ RU - PW) (2.4) 

G z + F z = m(W + PV - QU) (2.5) 

L » PI X - RJ XZ + QR(I Z -I y ) - PQJ XZ (2.6) 

M = QI y + RP(I X -I Z ) + (P 2 -R 2 )J x2 (2.7) 

N - ftl z - PJ XZ + PQ (I -I x ) + QRJ XZ (2.8) 

where , in body coordinates: 

U , V, W are the components of V 

P, Q, R are the components of w 


G x , G y , G z are the components of the gravity forces 

F x , F , F z are the components of the aerodynamic 

forces 

L, M f N are the components of the external moments 
I x , I y , I z are the mass momen s of inertia 

J X2 is the cross moment of inertia (J xy and J zy are 

zero because the X-Z plane is a plane of symmetry) . 

This body axis frame can be referenced to an inertial 
coordinate system by a translation of the center of mass 
and three rotations. Using the subscript "I" to identify 
the inertial coordinate system, "3" the instantaneous 


i 


I 

» 

i 



\ 

\ 

j 

i 

i 
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body axes coordinate system, and M l" and "2" intermediate 
coordinate systems, the transformation takes place in the 
following order: 

(1) starting with the "I" coordinate system, rotate 
about Zj through an angle y to establish a new 
coordinate system x^, y 1# z^. 

(2) Then rotate about y^ through an angle e to 
establish a new coordinate system X2, V2» Z2» 

(3) Finally, rotate about X2 through an angle $ to 
get to the aircraft fixed coordinate system 
*3' y 3' Z 3 • 

• ■ • 

From this description one can see that $, 0, y point 

along *3, y 2 , Zi respectively, whereas P, Q, R lie along 

x 3' y 3' z 3* They are related by the equations: 

• • 

P * ♦ - ¥ sin 0 

Q = 0 cos $ + $ cos 0 sin $ (2.9) 

R ■ V cos 0 cos <$ - (5 sin ♦ 

The external forces and moments come from two sources: 
9 rav ^ty and aerodynamic forces. Gravity, because it acts 
through the center of mass, produces no external moments. 

It does produce an external force that decomposes into: 

G x - - mg sin 0 

Gy ■ + mg cos 0 sin ♦ 

G z ■ mg cos 0 cos ♦ 

along the X3, yj, axes, respectively. 


( 2 . 10 ) 
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up to this point very few assumptions have been made 
and as a result the equations are all nonlinear. By 
assuming that the aircraft will only make small perturba- 
tions from some steady state condition (a reasonable 
assumption for a longitudinal tracking task) the equations 
can be greatly simplified. 


Let: 

U - U Q + u P 
V = V Q + v Q 
W * Wg + W R 


« p Q + p 0 - o 0 + e 

= Q 0 + q $ = $ 0 + 4> 

= R 0 + r V » Yq + * 


( 2 . 11 ) 


Where the "0" subscript refers to the trimmed value and 
the small letters are the perturbation variables. The 
following trimmed condition was chosen: 


U Q ■ 949 ft/sec 
■ 0 degree 
a Q - 3 degrees 

v 0 * W 0 “ P 0 " Q 0 “ R o “ *0 = *0 * 0 


( 2 . 12 ) 


The steady state control settings for this flight condition 
turn out to be: 

Throttle - 41,712 pounds thrust 

Elevator ■ - 6 degrees (2.13) 

and all other controls trimmed out to zero. The increment- 
al control variables are represented by 6 e , 6 t for 
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elevator and thrust. They represent changes from the 
above trimmed value. 

This same linearization technique can be applied to 
the aerodynamic forces and moments. This is done by 
separating them into a steady state term and linear terms 
due to the perturbation variables. For example: 


F v + 
x 0 


where a * 

o 



3F x 

3F x 

n 4 ± 

3F 

a + * 

au 

U + 3a 

0 

0 3 “ 

0 3 ® 

3F x 

X cx 



35 e 

0 c 




'o 



• 

0 = 

q • 




(2.14) 


The same is done for F^, F 


2 # L, M , and N . 

The partial derivatives are a function of the charac- 


teristics of the aircraft and the particular trimmed condi- 
tion at which they are evaluated. For the trimmed condi- 
tion selected, F x , F z , and M are independent of v , p, r, 

’!> , $, 6 a , 6 r while F y, L, N are independent of u, w, q, 9, 

6q, <5 t , $£. As a result, all the corresponding partial 
derivatives are zero. 

Combining equations 2.3 - 2.14, expanding the trigono- 
metric functions, and neglecting all second order terms 
(such as "qw H or "rv") , results in the following nine 
equations : 

p * ♦ - 0 sin 0 Q 
q - 9 

r * l cos 0* (2.15) 


l 


r 


i 


i 


mu = -mg0cos0 o +X u u+X a a+Xqq+X < .^6 e +X (5 ^6 t +X jS ^6 f 

m(v+O o r)-mg*cos0 o +Y v v+Y p p+Y r r+Y 6a « a +Y 5 ^5 r 

m(w-U o q)=-mgesin0 o +2 u u+Z a a+Z-a+Z q q+Z 6e 6 e +Z 5t 6 t +Z 5f 6 f 

P T x" f J xz' L v v+L P P+L r r +L 6 a 6 a +L 6 r 5 r 

qIy*M u u +M a a+M-a+ M qq+ M 5e 5e+M 5t 6t+M 6f 6f 

rl -pJ =N„v+N_p+N r+N. S-+N, 5 r 
z ^ xz v p r r °a a a r 

This aircraft is considered to be in straight and 
level unaccelerated flight and then to be disturbed by 
deflection of the elevator. This deflection applies a 
pitching moment causing a rotation which eventually causes 
a change in F x and F 2 , but does not cause a rolling or 
yawing moment or any change in F y ; thus p=R=*V=0 and the 
equations 2.4, 2.6, 2.8, may be eliminated, resulting 
in the following equations: 

q*0 

mu*-mg9+X u u+X a a+X . 6 e +X fi 

U 0 mi=n.U 0 9 + z u u+Z a a+ Zi i+Zj8+Z {e J e (2,1 

I y 5 “ M u u+M <i a+M i'’* M 9® +M «e s e 

where X f , Z, , Z. , M , end M are assumed zero, and 3 Q =0. 
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II. Elastic Mode Equations 

The most common method of accounting for the dynamic 
aeroelastic effects is to represent the deflection of the 
aircraft in terms of its normal modes of free vibration. 

The instantaneous deflection of any point on the structure, 
from the rigid steady-state condition as mtasured in the 
body axes frame, can be expressed as a vector sum of 
deflection components in the X, Y, and Z directions. 
Deflection in the Z direction and perpendicular to the 
XY plane is indicated by 5*(x, y, t) , where the planform 
has been idealized to a flat plate in the XY plane. 
Deflection in the Y direction and perpendicular to the 
XZ plane is £ y (x, z, t) , where the side view planform 
has been idealized to a plate in the XZ plane. Deflection 
in the X direction has a negligible effect on dynamic 

stability and is neglected. 

For longitudinal equations of motion we only need 

to consider the deflection in the Z direction which is 
£,(x, y, t) . The instenteneous deflections can be 
represented as a sum of an infinite number of orthogonal 
normal vibration modes as follows: 

QD 

l (x,y,t) - l ♦ i (x,y)n i (t) (2,1 

z i _ ^ 


i 


i 


i 
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where : 

$ (x,y) is the ith symmetric normalized mode shape 
• in the XY plane. 

n^(t) is the generalized displacement. 

Representation of vibration in terms of normal modes 
results in additional equations of motion of the form 
(reference 1) . 

m i C * n i + 2^ i w i n i +U) i n i 3 " (t) (2.18 

m ^ is generalized mass defined by 
mi - // m(x,y)$?(x,y)dxdy 
where : 

m(x,y) is the structural mass density per unit area 
in the respective XY idealized planform. 

(t) is the generalized force defined by 

Q n (t) - // f (x,y,t)$i (x,y)dxdy 
n i z 

where : 

f (x,y,t) is the time-dependent aerodynamic pressure 

z 

distribution acting in the Z direction. 

is an additional dependent variable to be intro- 
duced into the rigid-body equations as well as appearing 
in (2.18). Considering the longitudinal rigid-body 
equations (2.16), they will contain additional 


aerodynamic terms due to aeroelasticity as follows: 


L ni + ni l 

" hL. 3m_ • 

1 . snj. n i + 3n ± n i 

1®1 U A 1 J 


In order to make the elastic equations (2.18) com- 
patible with the rigid-body equations, the (t) must be 
expressed in terms of the dependent variables. Thus: 


+ I 3n. o k + l tt" \ 
k=l k k**l an k k 


( 2.20 


The numerous partial derivatives are obtained by 
utilizing the stability derivatives of the aircraft. 

Listed in Table 2.1 is the relationship between these 
partial derivatives and the stability derivatives while 
Table 2.2 gives the values of these derivatives for the 
B-l bomber at the Mach 0.85, sea level flight condition. 

Using Tables 2.1 and 2.2 to evaluate the coefficients 
in equations 2.16 and 2.18 results in the following 
equations , 
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u » -0 . 025U-25 . Oa-32.20+0.0001416^ 
a » -0.00065u-1.205o+0.943§-0.00905ni“0.00021n 1 

+0.00465n2 + 0* 00015 ^2"°* 28885 e 

6 = -0.0026u-7.649a-0. 481o-l. 5660-0. 1889ni“0.00754fii 
+O.210n 2 +0-00 704 n2“ 15 - 1936 e (2 2] 

n 1 +0.272n 1 +184.69n 1 — 735. 19o+2. 2640-135. 404fi+7.263ni 
-0.866^-15. 88 9n2 +0 - 915 ^2" 2228 * 4 5 e 

n 2 +0.424ri 2 +448.59n2” 764 . 70a+6 . 153a+50 . 1499+7 . 041ni 

- 0 . 1206 n 1 " 7 * 962 n 2 " 0 * 426 ^ 2 +614 * 96<S e 

Additional details of the derivation of perturbation 
equations for an aircraft simulation are presented in 


references 6, 7, 14. 


I 


I 




i 


I 
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Table 2.1 Force and Moment and Elastic Force Derivatives 
as a Function of Stability Derivatives 

a x - pUgS/2 

a 2 * a l c / 2U 0 


x u - ai c x /U 0 

U 

z u * a l c z/ D 0 

«u ■ a l cC m u ^ u 0 

*0 “ a l c x a 

z a * a l C z a 

M a ‘ a l cC m a 

Xa * a 2 c xJ, 

z i " a 2 C z- 

Ml = a,cC . 
° 2 “a 

X 8 - a 2 C x ‘ 

z e ‘ a 2 c zJ 

M 0 " a 2 cC mj 

x 6 e = a l C x Se 

z «e ‘ a l c z ie 

M Se “ a l c< V 

°e 

x St * 

2{ t = 0 

M 6t ' ° 

X n " a l C v 

^1 1 Xf U 

z n = a l C z 
n i Zn l 

M n " a l cC m n n l 
1 n l 

\ * a l C * • /U 0 

1 n i 

z n x * a l c z n -/ u o 

Mj. * a i cC % /u o 

% “ a l C *n 2 

Z = a.C 

n 2 * z n 2 

M » a.cC 

n 2 1 m n 2 

\ ’ aiCx 'i 2 /U ° 

Z • = a-iC • /U_ 

^2 ^ Zti 2 ® 

M-» ** a,cC . /U. 

*2 1 m n 2 o 


•Obtained from reference 3. 


i 
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Table 2.1, cont. 


° n l 

a,C 
1 "i 

S = 

a. C 

1 n 2 

a 

a 

a 

a 

Q 

n l: 

a 2 C n 
2 n le 

* 

a-,0 
2 n 2 • 

a 

a 

a 

a 

Q v 

a 2^"n 

*•« 

Q_ = 

n 2e 

a.,0 

2 n-j . 
^ 9 

Q »1 

3 a l c n 1 

°n 2 

=» a.C 
1 n 2 

”1 

n l 

n l 

n l 

V 

3 a i c ni . /u 0 


= a l C n 2 . /U 0 


*1 

n l 

n l 

N 

ln 2 

3 a l c ni 

n 2 

Q h 

n 2 

n 2 

3 a l Cn 2 

n 2 

° n « 2 

3 a i Cr ^ 2 /u o 

° n 2. 

n 2 

* a i%. /u o 
^n2 

Qn l«e 

= a, C 

1 1, 

<5e 

\ 

3 a l C n2. 

e 


♦Obtained from reference 3. 


i 
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Table 2.2 Stability Derivatives for B-l Bomber in 
Mach 0.85 Flight Condition 


C x - 
x u 

-0.08066 

c, » -1.9659 
z u 

C m u * 

-0.4546 

c ** 
x a 

-0.08500 

C_ - -3.9367 
a 

’ 

-1.41052 

c x . “ 

0 

C * -5 

z a 

c = 

m a 

-11.005 

C • * 
x 0 

0 

C-. - 17.8558 
z 9 

C • - 

m e 

-35.7556 

C *C 

0 

C z * -.9426 
5 e 

c “ 

-2.799 

CXn l 

= 0 

C_ =» 0.02922 

z 

n l 

C ®ni 

* -0.0348 

c • 

- 0 

C - 0.6592 

Z \ 

C • 
ra ni 

- -1.32169 

C xn 2 

“ 0 

C = -0.015 

Zn 2 

c 

m n 2 

= 0.0387 

C xn 2 

- 0 

C . =-0.4733 

Zn 2 

C mf. 
n 2 

= 1.233 
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Table 2.2, cont. 


C n ■ -0.06478 

U 1 

■••a 

C - 0.48975 

" 2 a 

C. - 0.02469 

n l. 

a 

C = 0.48779 

n 2 . 

C n ■ -1.47658 

1 • 

9 

C . = 3.97547 
n, • 

2 0 

C - 0.00064 

n l 

n l 

C * 0.00451 

n 2 

n i 

C„ » -0.07243 

C n = -0.0733 

2n l 

c n - -0.0014 

^ 2 

C n = -0.0051 

n 2 n 2 

C ni . - 0.0765 
n 2 

C n . = -0.2588 

n 2n 2 

c 0i Je “ -°- 19635 

C n = 0.3939 

n 25 e 

Note: (1) All coefficients are non-dimensional 

(2) Data obtained from 

reference 3 

m = 7085 slugs 

Uq = 949 ft/sec 

s ** 1946 ft 2 

p = 0.002378 slugs/ft 

c * 15.33 ft 

q « 1071.754 lbs/ft 2 

I y - 5.8916xl0 6 slug-ft 2 


2.3 Flight Director Equation 



The total pitch angle time history that the pilot 
feels and sees, either on the outside horizon or the 
attitude indicator display, is given by (reference 4 and 
5) . 
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2 

0 i (x p't) - 0(t) - J ♦i' (XpJnitt) (2.22) 

Figure 2.2 show3 the pitch angle response to 
command input due to rigid-body notion and that 
due to the two lowest frequency elastic mode 
responses. Where x^ indicates pilot fuselage station and 
4>^(Xp) the slope of the i th symmetric elastic mode at that 
station. 

The relative contribution of the elastic terms to 
the total pitch response increases as the natural fre- 
quencies of the modes decrease. Now the values of 
and $2 » which are obtained from B-l data (reference 3), 
are as follows: 

- 0.25 

= 0.29 

2 

So equation (2.22) can be rewritten as 

B t (x ,t) = 0(t) - O^Snj^ - 0.29n 2 (2.23) 

2.4 Pitch Angle Error Equation 

A pitch angle error results from the reference 
trajectory which is the commanded signal that the pilot 
is to follow. When the CRT is implemented, this pitch 
angle error would be a measured quantity. That is, the 
CRT would provide the aircraft with its present pitch 
angle and this would be compared with its desired pitch 
angle given by the reference trajectory. The difference 
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between the actual flight direction and the reference 
trajectory is the pitch angle error (see Figure 2.3). 

The equation is as follows: 

e 0 ** ©i - e c (2.24) 

The pilot's job is to maintain this pitch angle error as 
close to zero as possible. The electronic attitude dir- 
ector is the display presently on the simulator built by 
William Seitz when the simulator was located in the M.E. 
school at Purdue (reference 13) . A drawing of the dis- 
play is shown in Figure 2.3. The dashed medium length 
line is attached to the face of a cathode ray tube (CRT) 
and represents the wings of the aircraft. The single 
long line moves up and down due to aircraft motion and 
represents the horizon. When the aircraft is pitched up, 
the horizon line is below the fixed aircraft symbol. The 
two short lines on the EADI are often referred to as a 
flight director. They move up and down to give pitch 
commands. The pilot's job is to keep the aircraft symbol 
centered in the two short lines. This display can be 
configured in two ways. The first is called a control 
director. With the control director, movement of the 
horizon line is directly related to movement of the 
control surface (control stick) . When the pilot pulls the 
stick back (up elevator) the horizon line goes down. The 
displacement of the horizon line is related to control 
movements by a constant. As a result the two short lines 



Fixed Aircraft 


Command 


Horizon Lino 


(Aircraft la Pitched Up) 


(Rigid Body B c ) 


Flight Director 
(Command in g Pitch Down) 


Figure 2.3a Electronic Attitude-Director Indicator (EADI) 




Local H. 



Figure 2.3b The Airplane Attitude Correaponding to Above 

I 
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command elevator deflection. The other configuration 
ia a flight director. The horizon line and the two 
short lines are present in this mode also, they are just 
driven differently than the control director. A change 
in pitch angle (Q^) produces a vertical displacement of 
the horizon line. As a result, movement of the two short 
lines commands change in pitch angle. The electronics 
driving the display are presently configured as a flight 
director, rather that control director, we have used it 
as a flight director in this work. 

2.5 Summary Equations 

Equations 2.21 can be written in matrix form in 
Table 2.3 

Equations 2.21 are the longitudinal dynamics of the 
B-l bomber at mach 0.85, sea level flight condition 
without the structural mode control system operating. 

Equation 2.23 is the flight director equation and 
equation 2.24 is the pitch angle error equation. 

These equations are all linear and are simulated on 
an analog computer in the piloted simulation used in 


this study. 


Table 2.3 Equations of Motion in Matrix Form 


x = Ax + Bu 

Where: x = (u,a , 0 , 0 , nj # / n 2 » n 2 ) T 

u - (<5 e '<5 t )T 


-.025 

-25 

-32.2 

0 

0 

0 

0 

0 

-.00065 

-1.205 

0 

.943 

-.00905 

-.00021 

.00465 

.00015 

0 

0 

0 

1 

0 

0 

0 

0 

-.00228 

-7.0694 

0 

-2.0195 

-.1845 

-.00744 

.2078 

.00697 

0 

0 

0 

0 

0 

1 

0 

0 

-.00147 

-737.92 

0 

-133.268 

-177.447 

-1.1385 

-15.414 

.91534 

0 

0 

0 

0 

0 

0 

0 

1 

-.00399 

757.29 

0 

55.951 

6.985 

-.1219 

-465.523 

-.8491 


0 .000141 1 


-.288 0 

0 0 

-15.0541 0 

0 0 

-2229.054 0 

0 0 

613.183 0 


B * 
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CHAPTER 3 

DYNAMIC CHARACTERISTICS OF EIGHT CASES 
3.1 General Description 

In this study, the equation of motion as shown in 
equation (2.21) includes two elastic modes. The natural 
frequencies of the two elastic modes will be parametrically 
reduced. They are represented by eight cases which were 
each flown by several pilots. Dynamic characteristics 
for each of these eight cases will be specified by four 
modes: short period mode, phugoid mode, elastic mode 1, 

and elastic mode 2. The dynamics of these eight cases will 
be described in the following sections. 

3.2 n r .mlc Characteristics o f Case 1, Case 2, and Case _3 
In this section, only the natural frequency of elastic 
mode 1 <«!> was reduced and the natural frequency of 

elastic mode 2 <w 2 > wa * unchanged. 

The value of u x was varied from 13.59 rad/sec to 
4.24 rad/sec. Three cases were chosen in this range of 
Ul . The root- locus plot is shown in Figure 3.1. 

The value of ui is 13.59 rad/sec: it is the original 

value of wy. It was chosen as Case 1. 



Figure 3.1 The Locus of Roots of Matrix A, Case 1# Case 2 , 
Case 3, as a Function of Natural Frequency of 
Elastic Mode 1# 



I 
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As the value of was varied from 13.59 rad/sec 
to 9.21 rad/sec, all inodes (phugoid, short period, elastic 
mode 1, elastic mode 2) were stable. 

When the value of is less than 9.21 rad/sec, the 
phugoid mode became unstable. 

A value of at 9.17 rad/sec was chosen as Case 2. 

A value of at 6.16 rad/sec was chosen as Case 3, and 
the phugoid mode is unstable and non-oscillatory . This 
value of was the lowest value that could be chosen 
for the pilot to fly. When the value of was less than 
6.16 rad/sec, the pilot could not control the unstable 
phugoid mode due to the mode interaction with elastic 
mode 1. 

The dynamic characteristics of Case 1, Case 2, and 
Case 3 are shown in Tables 3.1, 3.2, and 3.3, respectively. 

3.3 Dynamic Characteristics of Case 4 

In this section, only the natural frequency of elastic 
mode 2 (d^, was reduced and the natural frequency of 
elastic mode 1 (u^) was unchanged. 

The value of w 2 was varied from 10.14 rad/sec to 
2. 83 rad/sec, and just one case was chosen in thi; range 
of <i)2 - The root-locus plot is shown in Figure 3.2. 

The value of uj 2 at 4.79 rad/sec was chosen as Case 4. 
This value of w 2 was as low as could be chosen as 
previously discussed. If the value of u 2 1* larger than 


I 
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Table 3.1 Dynamic Characteristics of Case 1 


Case 1 (original case) : 

* 13.59 rad/sec; 25 ]^ = 0.272 
w 2 = 21.18 rad/seci 2C 2 W 2 " 0,424 


Parameter 



Coupled ur arapec 
natural freq. 
(rad/sec) 


Coupled damping 
ratio __ 


Short Period 


-1.498 
+2 . 373 j 


2.806 


0.5339 


- 0.00139 
+ 0 . 00708 j 


.07081 


0.01973 



Coupled damped 
freq. (rad/s ec) 


Time to half or 
twice Amp. 

( sec ) 


Period 
(sec) 


2.373 


0.46 


0.0708 


494.7 












Elastic Mode 1 Elastic Mode^ 


- 0.6583 
+13 . 295 j 


- 0.^603 
+ 2 1 . 3 4 9 j 


13.312 


21.354 


0.0494 


0.0215 


13.296 


21.349 


1.49 


0.47 


0.29 


Table 3.2 Dynamic Characteristics of Case 2 


Case 2: 

W 1 = 9.17 rad/sec: 2 ^ 0 ^ = 0.183 
w 2 = 21.18 rad/sec; 2 q u - 0.424 

5 “> 


Parameter 

Roots 


Coupled undampe 
natural freq. 
(rad/sec) 

Coupled damping 
ratio 



Time to ha 
twice Amp. 
(sec) 


or 


Short Period 

Phugoid 

-1.3468 

♦2.1917j 

+3.4536xl0” 5 

+5. 7305x10-2 j 

2.5724 

0.057305 

0.5235 

-0.00060267 

2.19 

0.057304 

0.5124 

19979.149 

. 


Period 

(sec) 


2.8667 


109.6447 



KjJ 

ro 







Table 3.3 Dynamic Characteristics of Case 3 
Case 3: 

0 )^ = 6.16 rad/sec; 2C^w^ = 0.123 
u) 2 = 21.18 rad/sec; 2 ^ 2 U> 2 = 0,424 


Parameter 



Coupled undampec 
natural freq. 
(rad/sec) 


Coupled damping 
rdLic 



Time to ha 
twice Amp. 
(sec) 


Period 

(sec) 


Short Period 


-9.2283X10* 1 
+1 . 5093 j 


1.7691 


0.5217 


1.5093 


0.7476 


4.1631 


+9.0978x10 ; 
-7.6723X10" 2 



Elastic Mode 1 Elastic Mode 2 

-1.1728 " 4.5568x10 ^ 

+5 ! 7485 j + 2 1 . 3 5 2 j 


5.8669 


21.357 


0.1999 


0.021337 


5.7485 


21.3521 


0.5883 


1.5142 


1.0930 


0.2943 


U> 

CJ 




-.6 


Matrix A, Case 4, as a Function of 
Elastic Mode 2 , u> 2 
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4.79 rad/sec, the situations would be the same as Case 1 
and Case 2 . 

The dynamic characteristics of Case 4 are shown in 
Table 3.4. 

3.4 Dynamic Characteristics of Case 5 and Case 6 

Both natural frequencies of elastic mode 1 and 

elastic mode 2 (oj 2 ) were reduced. 

The values of and u 2 were varied from 13.26 rad/sec 
to 4.89 rad/sec. Case 5 and Case 6 were chosen in this 
range for and u>2. The root-locus plot is shown in 
Figure 3.3. 

When the values of a nd w 2 were rad/sec, the 

phugoid mode was unstable and oscillatory. This was 
chosen as Case 5. 

When the values of and w 2 were 6.93 rad/sec, the 
phugoid mode was unstable and non-oscillatory . This was 
chosen as Case 6. These values of and uj 2 are as low 
as could be chosen in Case 6, as previously discussed. 

The dynamic characteristics of Case 5 and Case 6 are 
shown in Tables 3.5 and 3.6, respectively. 

3.5 Dynamic Characteristics of Case 7 and Case 8 
The natural frequency of elastic mode 2 (^ 2 ) was 

varied from 10 rad/sec to 6.32 rad/sec while the natural 

frequency of elastic mode 1 (w^) was varied from 10 rad/sec 

to 12.61 rad/sec. The root-locus plot is shown in Fig- 
ure 3.4. 


Table 3.4 Dynamic Characteristics of Case 4 
Case 4: 


= 13.59 rad/sec? 2C 1 uj 1 = 0.272 
w 2 = 4.79 rad/ sec; 2c 2 w 2 « 0.096 


Parameter 

Short Period 

Phugoid 

Elastic Mode 1 

Elastic Mode 2 

Roots 

-1.0819 
+1 . 1 4 3 8 j 

+1.4654X10' 1 

-1.3167x10-1 

- 7.0119X10" 1 
+1 3 . 2 5 1 j 

-6.7877X10" 1 
+ 5 . 9 3 1 5 j 

Coupled undamped 
natural freq. 
(rad/sec) 

1.5745 


13.270 

5.9702 

Coupled damping 
ratio 

0.6872 


0.05284 

0.1137 

Coupled damped 
freq. (rad/sec) 

1.1438 


■ — 

13.251 

5.7544 

Time to half or 
twice Amp. 

(sec) 

0.6377 


9.8405 

r- — 

1.0165 

Period 

(sec) 

5.4931 


4.7415 

1.0919 













Table 3.5 Dynamic Characteristics of Case 5 


Case 5 : 


= 11.66 rad/sec; 25]^ 


0.233 


^2 = 11.66 rad/sec; 2£2 U 2 


0.233 


Parameter 

Short Period 

Phugoid 

Elastic Mode 1 

Elastic Mode 2 

Roots 

-1.4035 
+2. 1671 j 

+6 . 0288x10“** 
+5.3748xl0"2j 

- 9 . 1264x10“! 
+11 . 765 j 

- 1. 8770x10"* 1 
+11 . 573 j 

Coupled undamped 
natural freq. 
(rad/sec) 

2.5819 

0.053745 

11.801 

11.574 

Coupled damping 
ratio 

0.54359 

-0.00011217 

0.077338 

0.01621 

Coupled damped 
freq. (rad/sec) 

2.1671 

0.053745 

11.7657 

11.5725 

Time to half or 
twice Amp. 

(sec) 

0.4916 

114454.8772 

0.7560 

1 

3.6778 

Period 

(sec) 

2.8993 

116.9073 

0.5340 

0.5429 


































Table 3.6 Dynamic Characteristics of Case 6 


Case 6 : 


= 6.93 rad/sec; 2t^w^ =» 0.139 
u )2 * 6.93 rad/sec; 2c 2 w 2 = 0.139 


Parameter 

Short Period 

Phugoid 

Elastic Mode 1 

Elastic Mode 2 

Roots 

-9.0635X10" 1 

+1.7581x10'“ J 

-1.4072 

-5.2984xl0“ 2 


t9.7212xl0” 1 j 

-1.5307X10 -1 

+7 . 1942 j 

+ 6. 9716 j 

Coupled undamped 
natural freq. 
(rad/sec) 

1.3665 


7.3305 

6.9178 

Coupled damping 
ratio 

0.70279 


0.19196 

0.0075997 

Coupled damped 
freq. (rad/sec) 

0.9721 


7.1942 

6.9176 

Time to half or 
twice Amp. 

(sec) 

0.7185 


0.4903 

13.1246 

Period 

(sec) 

6.4634 


0.8734 

0.9083 




Figure 3.4 The Locus of Roots of Matrix A, Case 7, Case 8, as a 

Function of Natural Frequencies of Elastic Mode 1, top 
and Elastic Mode 2 , 
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When the values of and w 2 were 10.25 rad/sec and 
9.75 rad/sec, the phugoid mode and elastic mode 2 were 
unstable and oscillatory. This was chosen as Case 7. 

When values of and w 2 were 10.68 rad/sec and 9.27 

rad/sec, only the phugoid mode was unstable. This was 
chosen as Case 8. 

The reason that both Case 7 and Case 8 were chosen 
was that their dynamic characteristics were quite similar, 
yet elastic mode 2 in Case 7 was a little unstable and 
that in Case 8 was stable. What would the handling qualities 
and pilot ratings in these two cases be? 

The dynamic characteristics of Case 7 and Case 8 
are shown in Tables 3.7 and 3.8, respectively. 

3.6 Summary of Eight Cases 

All these eight cases included most of the situations 
in which the handling qualities and pilot rating would 

be affected differently by the two elastic modes included 
in the model. 

Each of the eight cases programmed in the analog 
computer driving the simulator displays was investigated 
to answer the questions posed. 

Case 1, Case 2, and Case 3 were chosen when only 
the natural frequency of elastic mode 1 was reduced. 

Case 4 was chosen when only the natural frequency 
of elastic mode 2 was reduced. 


Table 3.7 Dynamic Characteristics of Case 7 
Case 7: 

Wj » 10.25 rad/sec; 2^“^ =* 0.205 
u 2 * 9.75 rad/sec; 2C 2 t * , 2 “ 0,195 


Parameter 

Short Period 

Phugoid 

Elastic Mode 1 

Elastic Mode 2 

Roots 

-1.3207 
+1 . 9964 j 

+1.3608xl0" 3 
+2 . 8139xl0“ 2 j 

- 1.1554 
♦10.168j 

+ 4.2338X10" 3 

+ 9 . 8978 j 

Coupled undamped 
natural freq. 

/ rnA/&GC ) 

2.3937 

0.028172 

10.234 

9.8978 

Coupled damping 
ratio 

0.55174 

-0.048305 

0.1129 

- 0.00042775 

Coupled damped 
frpa. (rad /sec) 

1.9964 

0.02814 

10.1686 

9.8345 


0.5225 

507.0366 

0.5972 

162.9748 

Period 

(sec) 

3.1473 

223.2901 

0.6179 

0.6388 













Table 3.8 Dynamic Characteristics of Case 8 
Case 8: 

» 10.68 rad/sec; 26^ » 0.214 

w » 9.27 rad/sec; 26 2 w 2 “ °* 185 
2 


Parameter 



Short Period 


-1.3259 
+1. 9876j 


2.3893 


+1. 3833x10“-* 
+2 . 5533x10“^ j 


0.02557 


Coupled damping 
ratio 



Time to 
twice Amp. 
(sec) 


0.55493 


1.9877 


0.5204 


-0.054096 


0.02553 


498.8307 


Period 
(sec) 


3.1611 


246.0852 



Elastic Mode 1 

- 1.1404 
+10. 284 j 

10.347 

0.11021 

10.2839 

0.6051 


Elastic Mode 2 

- 5.1889xl0' 3 

+ 9 . 77 81 j 

9.7781 

0.00053066 

9.7781 

132.9775 


0.6109 


0.6426 
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Case 5 and Case 6 were chosen when both the natural 
frequencies of elastic mode 1 and elastic mode 2 were 
reduced. 

Case 7 was chosen when elastic mode 2 was unstable. 

Case 8 was chosen with its dynamic characteristics 
close to that of Case 7 but with elastic mode 2 stable. 

Table 3.9 shows the coupled undamped natural fre- 
quencies and damping ratios of the rigid and elastic modes 
for each of the eight cases. 






i 

1 £ 


Table 3.9 Natural Frequencies and Damping Ratios of Eight Cases 


Case 

# 

rad/ sec 

U)2 

rad/sec 

Ggp 

rad^Sec 

Cph ^Rh 

rad/sec 

*»le 

rad/sec 

*2e 

w 2e 

rad/sec 

1 

13.59 

21.18 

0.5339 

2.806 

0.0197 

0.0708 

0.0494 

13.312 

0.0215 

21.354 

2 

9.17 

21.18 

0.5235 

2.5724 

-0.00060267 

0.0573 

0.08769 

8.7891 

0.0213 

21.2'6 

3 

6.16 

21.18 

0.5217 

1.7691 

Real Roots 

+0.090978 

-0.076723 

0.1999 

5.8669 

0.0213 

21.357 

4 

13.59 

4.79 

0.6872 

1.5745 

Real Roots 

+0.14654 

-0.13167 

0.05284 

13.270 

0.1137 

5.9702 

5 

11.66 

11.66 

0.5436 

2.5819 

-0.0001122 

0.0537 

0.0773 

11.801 

0.0162 

11.574 

6 

6.93 

6.93 

0.7020 

1.3665 

Real Roots 

+0.17581 

-0.15307 

0.1919 

7.3305 

0.007599 

6.9178 

7 

10.25 

9.75 

0.5517 

2.3937 

-0.0483 

0.0282 

0.1129 

10.234 

-0.0004277 9.8978 

8 

10.68 

9.27 

0.5549 

2.3893 

-0.0541 

0.0256 

0.11021 

10.347 

0.0005306 9.7781 



46 


CHAPTER 4 
EXPERIMENT 


4.1 Overview 

A pilot in the loop simulation was performed to 
answer the following questions. Of the eight cases exam- 
ined, which one results in the worse pilot performance? 
What are the pilot comments on each case? How well can 
pilots perform in these cases compared to the original 
case (Case 1)? 

The eight cases previously discussed were used to 
answer these specific questions. They are. 

(1) Does Case 1 provide good handling qualities 
and pilot rating? 

(2) Are the handling qualities and pilot rating 
for Case 2 and Case 3 acceptable, in which only 
the natural frequency of elastic mode 1 is 
changed? The natural frequency of elastic 
mode 1 is changed to a greater extent in Case 

3 than that in Case 2. 

(3) Are the handling qualities and pilot rating for 
Case 4 acceptable, in which only the natural 
frequency of elastic mode 2 is changed? 
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(4) Are the handling qualities and pilot rating for 
Case 5 and Case 6 acceptable, in which both the 
natural frequencies of the elastic modes are 
changed? 

(5) Are the handling qualities and pilot rating for 
Case 7 and Case 8 acceptable? Elastic mode 2 

is slightly unstable in Case 7, but it is stable 
in Case 8. 

4 . 2 Apparatus 

A fixed base simulator was adapted to answer these 
questions. A block diagram of the apparatus involved in 
this test is shown in Figure 4.1. Two Applied Dynamics 
analog computers comprising 60 amplifiers were used to 
simulate the equations of motion, flight director 
equation, and command signal. The mathematical forms 
for these equations were presented in Chapter 2. A 
schematic diagram of this simulation is shown in Appen- 
dix. 

The cockpit mockup contained a control stick, a 
throttle, vertical velocity indicator, airspeed indicator, 
angle-of-attack indicator, altimeter, and the electronic 
attitude-director indicator (EADI). 

A control stick provides elevator input to the 
aircraft simulation. Springs produce a restoring force 
on the control stick to return it to a neutral position. 



Figure 4.1 Block Diagram of Apparatus 


00 
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The spring stiffness was suitable for pilots to fly 
(40 pounds for full elevator deflection) . The resulting 
gradient provides good centering of the stick but does 
not require excessive force. The stick neutral position 
can be adjusted by usinq the elevator trim crank on the 
cockpit floor. 

A throttle provides thrust input. Both of these 
inputs, elevator and throttle, are measured by the wiper 
of potentiometers mechanically connected to the movable 
controls. These measurements are then electrically 
connected to the input of operational amplifiers for use 
in the simulator. 

Initially the aircraft is flying at 949 ft/sec at 
constant altitude with 0 degree flight path angle and 
<iq = 3°, 0q » 0°, and up elevator 6 e = -6 degrees. 

The electronic attitude-direction indicator (EADI) 
vertical velocity, airspeed, altimeter, and angle-of- 
attack indicators are all visible. Description of the 
instruments is provided in reference 9. 

The vertical velocity and angle-of-attack indicators 
are driven by miniature direct current (DC) servo systems 
that convert an analog computer DC voltage into an 
angular rotation of the needle within the instrument. 

The EADI is produced by drawing special symbols on cathode 
ray tubes (CRT) . Details of the electronics are provided 
in reference 9. 


i 


i 
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4 . 3 Methodology 

Data was recorded on four pilots flying each of eight 
cases during two separate approaches (replications) . As 
a result, a total of 64 data runs were performed. This 
does not include all the time each pilot spent learning 
the simulator and practicing each case. This section 
describes the details of how the experiment was made 
including the cases arranged, instructions given, and 
data recorded. 

To minimize the possibility of learning trends 
affecting the results of this experiment, the eight cases 
were presented in a different order to each pilot. The • 
order they were presented in by means of random table / 
(reference 12) is shown in Table 4.1. 


Table 4.1 Order of Presentation of Cases 



Each of the four pilots was instrument rated, 
iheir experience varied from private, military, to 
commercial. The average total flight hours was more 
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than 1600 hour.. A detailed description of the simulator, 
each display, the task, and the purpose of the experiment 
was first distributed to each pilot. The task was 
described as “continuously maintain the tracki g 
of Pitch angle e 9 as close to zero as possible while 

using throttle to maintain the trim airspeed value. " 

. half an hour to one hour getting 

Each pilot spent from hair an n 

acquainted with the simulator. This involved learning the 
iocation and operation of each instrument and the dynamics 
and limitations of the aircraft. By this time each had 
a good understanding of the experiment and what was 

required. 

After each pilot had been acquainted with the instru- 
ments, a two hour session of simulated flight for the 
first eight cases was conducted. Fifteen minutes were 
used for each case. Ten minutes of the fifteen minutes 
were spent in practicing with the particular case being 
examined. During this ten minutes the pilot flew at 
least twice - one for no co-muand signal provided, another 
one for conmand signal provided. This command signal 
i. Shown in Figure 4.2. The deviation of pitch angle 

error was visible on the EADI; this acts as a means of 

. ...... iea-n*rc situation. When his 

determining the pilot s lea.n-ng 

• _ _ « ■- — • ' * zT ^ 11:2 • 

performance nai reaches a p-a ' 

was taken. After these eight cases were finished, they 
were repeated in the same order. Thus, sixteen 
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actual data runs were performed. 

Each run begins with the aircraft at 562.2 knots at 
6000 feet altitude. During the first 10 seconds, the 
pilot flies at a commanded 0.85 deg/sec pitch up. During 
the 10 second to 50 second period, the pilot makes any 
trim and power adjustments to keep the aircraft pitched 
up at 8.5 degrees. At 50 second, the pilot is given a 
command signal to start pitching down. After 20 seconds 
elapsed, the reference pitch angle is changed to “8.. '3 
degrees. From the 70 second to 110 second period, the 
pilot makes any trim and power adjustments to keep the 
aircraft pitched down at -8.5 degrees. At 110 seconds, 
the pilot flies a commanded 0.85 deg/sec pitch up. At 
120 seconds the experiment is completed. 

In computing work, the tracking error, which was 
recorded by FM recorder, was sampled ten times per second 
and converted to digital words. The root mean square 
(RMS) of the tracking error was then computed off-line. 
The analysis was done later off-line by using a mini- 
computer. Because each pilot by observing the EADI knew 
how much the aircraft symbol was off the command signal, 
he was always aware of how well he was doing. Also time 
histories of six pertinent parameters were recorded by 
three two-parameter strip chart recorders. 

Each pilot was given a questionnaire at the comple- 
tion of the experiment. It requested he list comments 
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on each case and assign a Cooper-Harper rating to each 
case. A tape recording was also made of verbal comments. 
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CHAPTER 5 
RESULTS 


5 . 1 Method of Analysis 

This chapter compares the performance in each of the 
eight cases, including an analysis of the mean tracking 
performance, a Cooper -Harper rating, and comments by the 
pilots. Root mean square (RMS) error of some quantity 
is frequently used to evaluate tracking performance. In 
this study, RMS pitch angle error is used. The magnitude 
of the RMS pitch angle error is an indication of how well 
the pilots followed the trajectory '’’he task as described 
to the pilots was to "maintain the pitch angle error as 
small as possible." As a result, average RMS error is 
the most appropriate measure of performance. 

One must be careful when comparing the mean perform- 
ance because of the experimental nature of obtaining the 
results and variability of the pilots. For example, 
suppose one obtained results where the mean performance 
of system A is one unit better than system B, but the 
standard deviation for each was 10 units. One could not 
say with much confidence that system A was better than 
system B. With such large variance in the results the 
mean of A could have accidentally been better than the 
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mean of B. The technique of Analysis of Variance is a 
method for determining the probability that the results 
were due to chance. 

Sixty-four observations of the RMS pitch angle error, 
Cjj , are shown in Table 5.1. The notation Cij is used 
to denote the RMS pitch angle error for each observation. 
Subscript i is the row number corresponding to pilot 
performance while subscript j is the column number corre- 
sponding to case number. Two sets of observations were 
made for each pilot, thus there are two rows corresponding 
to each pilot. Cj denotes the mean of RMS pitch angle 
error in j th case; C is the overall mean. Analysis of 
variance (reference 10 and 11) was performed on the 
tracking error of pitch angle and is shown in Table 5.2. 
Degree of freedom means freedom to vary. Degrees of 
freedom for total cases are the number of observations 
in total minus 1 or 63; degrees of freedom for "between" 
cases are the number of cases minus 1 or 7; degrees of 
freedom for "within" cases are the sum of the number of 
observations within each minus 1 or 56. Table 5.2 
indicates that the differences in performance among the 
cases were statistically significant at the 0.01 level. 

F ratio is a non-dimensional quantity that reflects the 
probability that the indicated results are due to chance. 
The F ratio is interpreted by use of the F table (refer- 
ence 11) . This table is entered with the number of 


Table 5.1 64 observations of RMS Pitch Angle Error 



Case 

Cl 

c 2 

C 3 

C 4 

— 

C 5 

C 6 

C 7 

C 8 

Pilot 


1 

2 

3 

4 

5 

6 

7 

8 


i 

0.0177 

0.0234 

0.0711 

0.0257 

0.0205 

0.1119 

0.0220 

0.0185 

P 1 

2 

0.0285 

0.0052 

0.1082 

0.0246 

0.0509 

0.1488 

0.00851 

0.0179 


3 

0.0053 

0.0038 

0.1061 

0.0309 

0.0274 

0.0651 

0.0141 

0.0052 

P 2 

4 

0.0113 

0.0256 

0.1077 

0.0291 

0.0226 

0.0796 

0.0369 

0.0339 


5 

0.0408 

0.0284 

0.0824 

0.0381 

0.0382 

0.1033 

0.0220 

0.0232 

P 3 

6 

0.0186 

0.0190 

0.0784 

0.0360 

0.0214 

0.1748 

0.0575 

0.0272 


7 

0.0192 

0.0240 

0.1121 

0.0395 

0.0180 

0.1791 

0.0261 

0.0199 

P 4 

8 

0.0191 

0.0166 

0.1258 

0.0396 

0.0114 

0.1943 

0.0190 

0.0160 


C^O.0201 C 2=0. 0183 C 3 =--0.09D0 C 4 =0.0331 ^=0.0263 C 6 =0.1321 C ? = 0 . 025 8 C 0 -O.O2O2 


C » .0481 
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Table 5.2 Analysis of Variance 


Source of 
Variation 

Degree of 
Freedom 


RMS Pitch 
Angle Error 



Sum of 
square 

Mean 

square F-ratio 

Between cases 

7 

0.10637 

0.01518 35.201** 

Within cases 

56 

0.02425 

0.00043 

Total 

63 

0.13042 


where: ‘‘significant at .01 

level 



8 8 

"Between" sum of square: J I ( C 4 ~ c ) 

i=l j-1 1 


8 8 2 

"Within" sum of square: J J ( c ij “ c *) 

i=l j=l J 

8 8 2 

"Total" sum of square: J I 

i«l j=l 13 

"mean squares" are obtained by dividing each of the 
sum of squares by its respective number of degrees 
of freedom 

mean square for "between" cases 
"F-ratio" ■ " 


mean square for "within" cases 


degrees of freedom for the greater mean square across the 
top and with the number of degrees of freedom in the 
lesser mean square on the left-hand side. For this 
problem, we go over to 7 and down to 56. In that loca- 
tion we observe that the value of F needed for signifi- 
cance at 1 percent point is 2.98. Since our obtained F 
is greater than this, this means that the chance that 
these means of cases are significantly different is 99 
percent. 

It should be emphasized that statistical significance 
and engineering significance are two different concepts. 
Statistical significance refers to the probability that 
the results are due to chance whereas engineering signifi- 
cance refers to the benefits due to increased performance. 
Statistical significance can be well quantified whereas 
engineering benefit is frequently very subjective. 

Unfortunately, Analysis of Variance tells only if 
the differences among C^, C 2 , C 3 , C 4 , C 5 , Cg, C 7 , and 
C 8 are statistically significant. What is of interest 
is a comparison among the cases. This can be done by 
using the "Newman-Keuls Test" (reference 12) . This test 
determines how large a difference is required between two 
means for this difference to be statistically significant. 
The difference reflects the magnitude of the means and 
variances and the number of pilots and replications 
performed. That is, as the means or variances increase. 
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so must the difference between the means to maintain the 
same level of significance. 

5.2 Comparison of Cases 

Table 5.3 gives the means and standard deviations for 
each of the different cases among all pilots. For example, 
the average RMS pitch angle error for the four pilots 
over two replications C^ is .0201 radians. These results 
are plotted in Figure 5.1. The diamond represents the 
mean RMS pitch angle error and the dashed lines extend to 
the one sigma deviations. 

Table 5.4 is the result of the "Newman-Keuls Test" 
discussed previously. It presents the required differ- 
ences between means for various levels of significance. 
Comparing the difference seen in Table 5.3 with those 
required as presented in Table 5.4 indicates that Cases 
Cg and C 3 were each significantly worse than C 4 , C 5 , 

C-j, Cg, C±, C 2 . In between Cg and C3, Cg was signifi- 
cantly worse than C3. While C4 appears to be worse than 
C5, Cy, Cg, C^, C 2 * the difference was not statistically 
significant. The above indicates that the pilots do 
equally well for cases C 4 , C 5 , C7, Cg, C^, and C 2 in 
controlling the pitch angle error. In eight cases, Cg 
is the worst, C3 is the second worst. 
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Table 5.3 Summary of Results 

RMS pitch 
Angle Error 

Radian (Degree) 




Mean 


Std. 

Dev 

Case 

C 1 

.0201 

(1.1516) 

.0107 

(.6131) 

Case 

C 2 

.0183 

(1.0485) 

.0093 

(.5329) 

Case 

C 3 

.0990 

(5.6723) 

.0192 

(1.1001) 

Case 

c 4 

.0331 

(1.8965) 

.0060 

(.3438) 

Case 

c 5 

.0263 

(1.5069) 

.0126 

(.7219) 

Case 


.1321 

(7.5688) 

.0488 

(2.7960) 

Case 

c 7 

.0258 

(1.4782) 

.0153 

(.8766) 

Case 

c 8 

.0202 

(1.1574) 

.0084 

(.4813) 
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Table 5.4 Required Difference 

Significance 

Level of 
Significance 

Cases 

0.01 level 


in Means for Statistical 


RMS 

Pitch Angle Error 
radian (degree) 
0.0311 (1.78) 


The pilots were asked to assign to each case a Cooper- 
Harper rating for the tracking task and the aircraft as 
Simulated. The means and standard deviations of these 
ratings are shown in Figure 5 . 2 . c X . C 2 , c 5 , and C a were 
good; c 7 required minimal pilot compensation; C 4 required 
minimal to moderate pilot compensation; c 3 required con- 
siderable to maximum compensation; c 6 required extensive 
to unacceptable compensation. This result is consistent 
With the tracking performance just presented. That is, 

C 6 is worse than C 3 which is much worse than c 4 , C ? , Cg , 

C 5 , C 2 , and C 3 . Here, c 4 is worse than C 7 which is 
worse than C 9 , c s . C 2 and C 3 . c x . c 2 . c,. and C 8 show 
little difference. 


Another important criterion for determining the 
"worst" case is the pilot's preferences and criticisms 
about the cases. In response to the request, "compare 
the difficulty of flying case 1 with the difficulty of 
the other seven cases,” the pilots' answers 


around: 


centered 


PILOT RATING 





PILOT 

PILOT 

COMPENSATION 

RATING 

REQUIRED 

8 

UNACCEPTABLE 

7 

MAXIMUM 

6 

EXTENSIVE 

5 

CONSIDERABLE 

4 

MODERATE 

3 

minimal 

a 

UNNEEDED 


Figure 5.2 Cocper-Harper Ratings 


o> 



65 


Ci » easy 

C2 = about same as C]_ 

C 3 * difficult 

C 4 » slightly difficult 

C5 = about same as 

C 6 - more difficult 

C 7 - slightly difficult 

C 3 * about same as 

in addition, the pilots gave the following comments on 
the eight cases: 
c l : 

( 1 ) Nothing objectionable on the pitch altitude 
tracking at all. Very nice case. 

( 2 ) It was very easy to hold it precisely. No 
control effort at all. 

(3) I felt this was the easiest case to con- 
trol. 

(4) I find no problem with the dynamic charac- 
teristics of this simulator. 

c 2 : 

(1) Fairly easy one to fly-it was not diffi- 
cult to hold the command symbol right with 
the airplane symbol. There was very little 
oscillation in it: very little control 
motion is required to follow the command. 
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(2) Pitch-wise a little more difficult than 

situation, but still able to maintain 
pretty much within the limits. 

(3) The slight continuous oscillation plus the 
slight lagging control response would be 
acceptable in a long flight in the event 
that it did not get worse. 

(1) It did appear to be unstable in rigid- 
body dynamics. It was fairly difficult 

to hold the airplane symbol in the command 
box. It took fairly strong amounts of con- 
centration to keep the symbol in the com- 
mand box and if you got away from it a 
little way, it was very difficult to get 
it back. 

(2) Unable to maintain within the pitch limits 
and had a problem of pilot-induced oscilla- 
tions due to large elastic amplitude 
oscillations . 

(3) This case was very objectionable due to 
the extreme lag in the change of pitch of 
the aircraft following a control move- 
ment. 


x 
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(1) Not quite as easy as some of the others, 
but it was not difficult at all either. 

Did seem to be a little unstable in rigid- 
body dynamic motion. It took a little 
more forward pressure to do it and stick 
has a very stiff spring in the forward 
motion. 

(2) The pitch control was a little bit annoy- 
ing. There did not seem to be quick 
enough response with the cyclic stick. 

c 5 - 

(1) That was about exactly the same difficulty 
as Case 2? that means its very easy to 

fly and was exactly the same in all charac- 
teristics as the Case 2. 

(2) Handling characteristics seem to be pretty 
good there - not too much excursion on 
pitch. 

(3) Has very good pitch control. 

c 6 5 

(1) That wab tremendous amount of oscillation 
in it. The control motions were very 
exaggerated and it took large displacement 
on the control to try to follow up. 


i 
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(2) Difficulty maintaining within the limits 
pitch wise, probably out on pitch. 

(3) The cyclic response is not real good. 

Abrupt immediate large changes produces 
disastrous effects. 

(4) With the severity of these oscillations, 
caused in this case, it would be virtually 
uncontrollable . 

(1) It was not terribly difficult to track 
the command box. There was noticeable 
oscillation, due to the elasticity; but 
again, it was not too difficult to ignore 
that and to fly simpiy the rigid portion 
of the pitch profile. 

(2) That was easier to control. The oscilla- 
tion due to elasticity was of main annoy- 
ance . 

(3) In this case, the aircraft would be quite 
controllable without an augmentation. 

(1) I believe that this case is the same as 
Case 2 and Case 5 - not difficult at all 
to fly the command profile. This one I 
can hold almost at the center all the time 
precisely. There was a little more 
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oscillation involved due to elasticity 
apparently, but it was high enough 
frequency that it was easier to ignore 
that and simply to fly the rigid-body of 
the profile that is coming from the 
command . 

All three performance measures, tracking error, 
Cooper-Harper Rating, and pilots' comments agree that 
C 6 was generally worse than C 3 which was worse than C 4 , 

C 7 1 Cg, Cg, C2» and C3. 

Table 5.5 shows the summary for each of eight cases 
in three performance measures. 

5 . 3 Time Histories 

A sample time history of six pertinent variables 
( 9 i » e Q/ ^ 2 ^ 2 '' 9' <5 e ) examined for each of the 

eight cases to see the how and why of the differences in 
pitch angle error. All the sample time history plots did 
not come from one pilot. 

Figure 5.3 is a sample time history plot for Case 1. 
During the 120 seconds flying time, is almost the same 
as 9, due to the fact that j. ' and n 2 ^ 2 ' are small and 
do not contribute much to 8 ^. This case approaches that 
of a rigid-body, eg was easy to zero out; that is, the 
dynamic characteristics are good. 
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Table 5.5 Summary of Tracking Error 


f 


* 



No . of 
Case 

Tracking 

Error 

Pilot 

Rating 


Radian 

Degree 

C 1 

.0201 

1.1516 

1.6 


C 2 

.0183 

1.0485 

2.0 


C 3 . 

. 0990 

5.6723 

5.9 


C 4 

. 0331 

1.8965 

3.1 


C 5 

. 0263 

1.5069 

2.0 


C 6 

. 1321 

7.5688 

6.7 


C 7 

. 0258 

1.4782 

2.3 


C 3 

. 0202 

1.1574 

1.9 



Cooper-Rating, and Pilots' Comments 

Pilots' Comments 
Good 

Very little oscillation 

Large amplitude oscillation, difficult 
control 

Not very easy to control 
A little more oscillation than C 2 
Tremendous amount of oscillation 
An annoying oscillation 
A little more oscillation than C 2 
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0.0116 



~ 0.0058 
■o 

2, o.oooo 

- e ?- 0.0058 
^ - 0.0116 


m 


Samp! 







Future work should use a full state feedback 
control law to place the roots of the characteristic 
equation at precise values for each case. The rigid- 
body dynamics could then be maintained at their nominal 
values and the elastic mode coupled frequencies placed 
where desired. This would ensure that the pilot ratings 
would be based on the relative amplitudes of rigid and 
elastic pitch angle responses as presented on the EADI 
and not on poor rigid-body dynamics. 
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APPENDIX 


ANALOG SIMULATION 


t 


The equations presented in Table 2.3 and equation 
(2.23) and (2.24) were normalized by the scale factors 
given in Table A.l and implemented on two Applied 
Dynamics analog computers. Patching diagrams of these 
equations are shown in Figures A.l - A. 8. Included are 
the longitudinal equations of motion, flight director 
equation, pitch command input, and control inputs. 

These differential equations were used to obtain the 
tracking error of pitch angle, angle-of-attack, airspeed, 
vertical speed, and altitude. Table A. 2 gives the poten- 
tiometer settings used for this simulation, and Table 2.3 
gives the potentiometer settings for each of the eight 


cases. 
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i 
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Figure 5.4 is a sample time history plot for Case 2. ! 

In this case is slightly greater than that in Case 

1. However, its contribution to Q i is still small. e 9 
can still be zeroed out easily. 

Figure 5.5 is a sample time history for Case 3. In ] 

this case is much greater than that in Case 2, 

causing greater contribution to 0 i . Significant changes 
in 9^ are observed. The pilot has difficulty zeroing out ! 

e^; eg is large. The difficulty in operating this case 

can be seen from the great movement of the control stick ! 

(as shown by 5 e plot) . Dynamic characteristics are poor ' 

in this case. 

Figure 5.6 is a sample time history plot for Case 4. ■ 

In this case the major contribution to 0^ is n 2 ^ 2 * " How_ 
ever, the effect of n 2 ^ 2 ' l ess than that of in 

Case 3; 0j_ in this case is smaller than 6j_ in Case 3. 

Even though the natural frequency of elastic mode 2 is f 

reduced to a greater extent than that cf elastic mode 1 ; 

in Case 3, the effect on 8^ is less. e 9 is slightly ; 

greater than that in Case 2. The pilot has to work harder 
to control the pitch angle (as shown by 6 e plot) . 

Figure 5.7 is a sample time history plot for Case 5. ■ 

Both an ^ n 2 ^ 2 ' contribute to ®i* However, their - 

magnitudes are small. The pilot still can zero it out; 

e 9 is not quite as large, eg has an oscillation of small | 

magnitude due to the interaction between and n 2<t 2 * * 

r 1 i m 5 j 
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This oscillation causes annoyance to the pilot but is of 
no major concern due to its small magnitude. 

Figures 5.8 is a sample time history plot for Case 6. 
In this case both and n 2 * are rauch 9 reater than 

that in Case 5, causing greater contribution to 9^ 6^^ 

is significantly increased. The pilot can not zero it 
out; eg is large. The difficulty in operating this case 
can be seen from the greater movement of the control 
stick (as shown by <5 e plot) . Dynamic characteristics 
are the worst among the eight cases examined. 

Figure 5.9 and Figure 5.10 are sample time histories 
for Case 7 and Case 8. anc ^ ^2^2* ^ ase ^ have 

greater contribution to 0j_ than that in Case 8. The 
phenomena of oscillation in ^ 2 ^ 2 ' P resen *- 

in both cases. Elastic mode 2 is slightly unstable in 
Case 7, while it is stable in Case 8. However, in both 
cases, the pilot car. operate the aircraft with ease, 
e 6 is small in both cases. 
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Figure 5.9A Sample Tine History - Case ? 










CHAPTER 6 


SUMMARY, CONCLUSIONS, AND 
RECOMMENDATION 


6 . 1 Summary 

The longitudinal tracking task was evaluated using 
a fixed base simulation of the aircraft equations of 
motion. Pilots' comments on the task revealed the effect 
of low frequency dynamic aeroelasticity on the handling 
qualities and the pilot rating. 

In Case 3 and Case 6 there is a large over shoot in 
the pitch angle error. These two cases received worse 
Cooper-Harper ratings indicating significant difficulty 
due to mode interaction. In Case 5, Case 7, and Case 8 
there is not too much error in the pitch angle and the 
Cooper-Harper ratings are fair, but the pilots' comments 
are that it was annoying due to the oscillation from 
aeroelasticity. For the rest of the cases, the aero- 
elasticity has little effect on handling qualities and 
pilot ratings because of the minimal interaction with 
the rigid-body dynamics. 


6.2 Conclusions 


All conclusions are restricted to a constant altitude 
trim flight condition, which is the only flight condition 
which was investigated. 

The most obvious conclusion is that the natural 
frequencies of elastic modes that were investigated repre- 
sent good to very bad pilot ratings, primarily because of 
adverse aeroelastic mode interaction effects with rigid- 
body dynamics. 

Elastic mode 1 affects handling qualities and pilot 
rating more than elastic mode 2 when the natural frequen- 
cies of the two elastic modes are reduced to the same 
numerical value. 

From the point of view of handling qualities of 
rigid-body dynamics, Case 3 seems to have better handling 
qualities than Case 4. However, it is not true. So we 

I I 

have to look at and ti 2 $2 in pilot time history 

for Case 3 and Case 4, respectively. The 0^ was affected 
more in Case 3 than in Case 4 since contributes 

more to the value of 0^ in Case 3 than n 2 *2 contributes 
to the ^alue of 0j_ in Case 4. 

As seen in Case 7 and Case 3, the pilot rating is 
the same as the original case (Case 1) . From the pilot 
time history and comments we know the small elastic 
oscillation can be visually separated from rigid-body 


dynamics by the pilot. The elasticity is only slightly 
annoying to the pilot since the elastic modes do not 
significantly affect the rigid-body dynamics. 

In Case 6 both natural frequencies of the two elastic 
modes were reduced to a still greater extent. This is 
the worst investigated and resulted in the poorest pilot 
ratings. 

The results of the study indicate that handling 
qualities and pilot rating are functions of the natural 
frequencies of the elastic modes. The lower the value 
of natural frequency of an elastic mode, the worse the 
handling qualities and pilot ratings. This was due to 
adverse coupling effects on the phugoid dynamics and a 
larger contribution from the elastic mode amplitudes 
to the total pitch angle 9^. 

6.3 A Recommendation 

In this study, the natural frequencies of the elastic 
modf.s were lowered from their nominal values to values 
close to the nominal short period frequency. However, 
this resulted in mode interaction which lowered the 
short period frequency and caused the phugoid mode to 
split into positive and negative real roots. Much of 
the pilots’ difficulty in tracking on the worst cases 
was thus due to the positive phugoid root. 


Table A.l Normalization of Variables 


Variable 

u 

a 

0 

• 

9 

8i 

n l 

h 

^2 

*2 

h 


. ^realization 
10u volts * 
400 ft/sec 

0.2 rad 

0.6 rad 

0.2 rad 

0.6 rad 

10 

50 

10 

50 

6000 ft 

1.0 rad 

50.000 lbs 

















Figure A. 6 Pitch Command Input 
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Table A. 2 

Potentiometer Settings 


2A1 

.144 

2B1 

.123 

2A2 

.029 

2B2 

.500 

2A3 

.130 

2B3 

.032 

2A4 

.121 

2B4 

.310 

2A5 

.943 

2B5 

.224 

2A6 

.453 

2B6 

.140 

2A7 

.053 

2B7 

.122 

2A8 

.233 

2B8 

.913 

2A9 

.038 

2B9 

.849 

2A10 

.013 

2B10 

.048 

2C1 

.446 

2D1 

.753 

2C2 

.500 

2D 2 

.333 

2C3 

.012 

2D3 

.456 

2C4 

.295 

2D4 

.707 

2C5 

.533 

2D5 

.202 

2C6 

.355 

2D6 

.923 

2C7 

.114 

2D7 

.186 

2C8 

.308 

2D8 

.104 

2C9 

.915 

2D9 

.174 

2C10 

.025 



1A1 

.625 

1F1 

.052 

1A2 

.725 

1F2 

.052 
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Table 2.3 

Potentiometer 

Settings for 

Eight Cases 

Case 1: 


Case 5: 


2C6 

.355 

2C6 

.258 

2C7 

.114 

2C7 

.110 

2B8 

.913 

2B8 

.288 

2B9 

.849 

2B9 

.658 

Case 2: 


Case 6: 


2C6 

.154 

2C6 

.082 

2C7 

.105 

2C7 

.101 

2B8 

.913 

2C8 

.112 

2B9 

.849 

2C9 

.564 

Case 3: 


Case 7 : 


2C6 

.062 

2C6 

.196 

2C7 

.099 

2C7 

.107 

2B8 

.913 

2B8 

.206 

2B9 

.849 

2B9 

.620 

Case 4: 


Case 8 s 


2C6 

.355 

2C6 

.214 

2C7 

.114 

2C7 

.108 

2B8 

.062 

2B8 

.188 

2B9 

.521 

2B9 

.611 


